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Abstract 

Our numerical simulations show that the reconnection of magnetic field becomes fast in the presence of weak tur- 
bulence in the way consistent with the Lazarian & Vishniac (1999) model of fast reconnection. We trace particles 
within our numerical simulations and show that the particles can be efficiently accelerated via the first order Fermi 
acceleration. We discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays 
measured by Voyagers. 
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1. Magnetic Reconnection in CoIIisionless and Col- 
lisional Fluids 

A magnetic field embedded in a perfectly conduct- 
ing fluid preserves its topology for all time (Parker 79). 
Although ionized astrophysical objects, like stars and 
galactic disks, are almost perfectly conducting, they 
show indications of changes in topology, "magnetic re- 
connection", on dynamical time scales (Parker 1970, 
Lovelace 1976, Priest & Forbes 2002). Reconnection 
can be observed directly in the solar corona ( Innes et al 
1997, Yokoyama & Shibata 1995, Masuda et al. 1994), 
but can also be inferred from the existence of large scale 
dynamo activity inside stellar interiors (Parker 1993, 
Ossendrijver 2003). Solar flares (Sturrock 1966) and 
y-ray busts (Fox et al. 2005, Galama et al. 1998) are 
usually associated with magnetic reconnection. Previ- 
ous work has concentrated on showing how reconnec- 
tion can be rapid in plasmas with very small collisional 
rates (Shay et al. 1998, Drake 2001, Drake et al. 2006, 
Daughton et al. 2006), which substantially constrains 
astrophysical applications of the corresponding recon- 
nection models. 

We note that if magnetic reconnection is slow in some 
astrophysical environments, this automatically means 
that the results of present day numerical simulations in 
which the reconnection is enevitably fast due to numer- 
ical diffusivity do not correctly represent magnetic field 
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dynamics in these environments. If, for instance, the 
reconnection were slow in collisional media this would 
entail the conclusion that the entire crop of interstellar, 
protostellar and stellar MHD calculations would be as- 
trophysically irrelevant. 

Here we present numerical evidence, based on three 
dimensional simulations, that reconnection in a turbu- 
lent fluid occurs at a speed comparable to the rms ve- 
locity of the turbulence, regardless of either the value of 
the resistivity or degree of coUisionality. In particular, 
this is true for turbulent pressures much weaker than the 
magnetic field pressure so that the magnetic field lines 
are only slightly bent by the turbulence. These results 
are consistent with the proposal by Lazarian & Vishniac 
(1999, henceforth LV99) that reconnection is controlled 
by the stochastic diffusion of magnetic field lines, which 
produces a broad outflow of plasma from the reconnec- 
tion zone. This work implies that reconnection in a tur- 
bulent fluid typically takes place in approximately a sin- 
gle eddy turnover time, with broad implications for dy- 
namo activity (Parker 1970, 1993, Stix 2000) and par- 
ticle acceleration throughout the universe (de Gouveia 
dal Pino & Lazarian 2003, 2005, Lazarian 2005, Drake 
et al. 2006). 

Astrophysical plasmas are often highly ionized and 
highly magnetized (Parker 1970). The evolution of the 
magnetic field in a highly conducting fluid can be de- 
scribed by a simple version of the induction equation 
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where B is the magnetic field, v is the velocity field, 
and T] is the resistivity coefiicient. Under most circum- 
stances this is adequate for discussing the evolution of 
magnetic field in an astrophysical plasma. When the 
dissipative term on the right hand side is small, as is 
implied by simple dimensional estimates, the magnetic 
flux through any fluid element is constant in time and 
the field topology is an invariant of motion. On the 
other hand, reconnection is observed in the solar corona 
and chromosphere (Innes et al. 1997, Yokoyama & 
Shibata 1995, Masuda et al. 1994, CiaraveUa & Ray- 
mond 2008), its presence is required to explain dynamo 
action in stars and galactic disks (Parker 1970, 1993), 
and the violent relaxation of magnetic fields following a 
change in topology is a prime candidate for the acceler- 
ation of high energy particles (de Gouveia Dal Pino & 
Lazarian 2003, henceforth GL03, 2005, Lazarian 2005, 
Drake et al. 2006, Lazarian & Opher 2009, Drake et al. 
2010) in the universe. Quantitative general estimates for 
the speed of reconnection start with two adjacent vol- 
umes with different large scale magnetic field directions 
(Sweet 1958, Parker 1957). 

The speed of reconnection, i.e. the speed at which 
inflowing magnetic field is annihilated by ohmic dissi- 
pation, is roughly 77/ A, where A is the width of the tran- 
sition zone (see Figure 1). Since the entrained plasma 
follows the local field lines, and exits through the edges 
of the current sheet at roughly the Alfven speed, Va, 
the resulting reconnection speed is a tiny fraction of the 
Alfven speed, Va = BjiAnpyi^ where L is the length of 
the current sheet. When the current sheet is long and the 
reconnection speed is slow this is referred to as Sweet- 
Parker reconnection. Observations require a speed close 
to Va, so this expression implies that L ~ A, i.e. that the 
magnetic field lines reconnect in an "X point". 

The first model with a stable X point was proposed 
by Petschek (1964). In this case the reconnection speed 
may have little or no dependence on the resistivity. The 
X point configuration is known to be unstable to col- 
lapse into a sheet in the MHD regime (see Biskamp 
1996), but in a collisionless plasma it can be maintained 
through coupling to a dispersive plasma mode (Stur- 
rock 1966). This leads to fast reconnection, but with 
important limitations. This process has a limited as- 
trophysical applicability as it cannot be important for 
most phases of the interstellar medium (see Draine & 
Lazarian 1998 for a list of the idealized phases), not to 
speak about dense plasmas, such as stellar interiors and 
the denser parts of accretion disks. In addition, it can 
only work if the magnetic fields are not wound around 
each other, producing a saddle shaped current sheet. In 
that case the energy required to open up an X point 



is prohibitive. The saddle current sheet is generic for 
not parallel flux tubes trying to pass through each other. 
If such a passage is seriously constrained, the magne- 
tized highly conducting astrophysical fluids should be- 
have more like Jello rather than normal fluids. 

Finally, the traditional reconnection setup does not 
include ubiquitous astrophysical turbulenc^j (see Arm- 
strong, Rickett & Spangler 1994, Elmegreen & Scalo 
2004, McKee & Ostriker 2007, Haverkorn, Lazarian 
2009, Chepurnov & Lazarian 2010). Fortunately, this 
approach provides another way of accelerating recon- 
nection. Indeed, an alternative approach is to consider 
ways to decouple the width of the plasma outflow region 
from A. The plasma is constrained to move along mag- 
netic field lines, but not necessarily in the direction of 
the mean magnetic field. In a turbulent medium the two 
are decoupled, and fluid elements that have some small 
initial separation will be separated by a large eddy scale 
or more after moving the length of the current sheet. As 
long as this separation is larger than the width of the 
current sheet, the result will not depend on 77. 

LV99 we introduced a model that included the effects 
of magnetic field line wandering (see Figure 1). The 
model relies on the nature of three-dimentional mag- 
netic field wandering in turbulence. This nature is differ- 
ent in three and two dimensions, which provides the ma- 
jor difference between the LV99 model and the earlier 
attempts to solve the problem of magnetic reconnection 
appealing to turbulence (Matthaeus & Lamkin 1985). 
The eff'ects of compressibility and heating which were 
thought to be important in the earlier studies (Matthaeus 
& Lamkin 1985, 1986) do not play the role for the LV99 
model either. The model is applicable to any weakly 
perturbed magnetized fluid, irrespectively, of the degree 
of plasma being collisional or coUisionless (cf. Shay et 
al. 1998). 

Two eff'ects are the most important for understanding 
of the nature of reconnection in LV99. First of all, in 
three dimensions bundles of magnetic field lines can en- 
ter the reconnection region and reconnection there inde- 
pendently (see Figure 1), which is in contrast to two di- 
mensional picture where in Sweet-Parker reconnection 
the process is artificially constrained. Then, the nature 
of magnetic field stochasticity and therefore magnetic 
field wandering (which determines the outflow thick- 
ness, as illustrated in Figure 1) is very different in 2D 



' The set ups where instabihties play important role include Simizu 
et al. (2009a,b). For sufficiently large resolution of simulations those 
set-ups are expected to demonstrate turbulence. Turbulence initia- 
tion is also expected in the presence of plasmoid ejection (Shibata & 
Tanuma 2001). Numerical viscosity constrains our ability to sustain 
turbulence via reconnection, however. 
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Figure 1 : Upper plot: Sweet-Parker model of reconnection. The out- 
flow is limited by a thin slot A, which is determined by Ohmic dif- 
fusivity. The other scale is an astrophysical scale L 2> A. Middle 
plot: Reconnection of weakly stochastic magnetic field according to 
LV99. The model that accounts for the stochasticity of magnetic field 
lines. The outflow is limited by the diffusion of magnetic field lines, 
which depends on field Une stochasticity. Low plot: An individual 
small scale reconnection region. The reconnection over small patches 
of magnetic field detemiines the local reconnection rate. The global 
reconnection rate is substantially larger as many independent patches 
come together From Lazarian et al. 2004. 



and the real 3D world (LV99). In other words, removing 
artificial constraints on the dimensionality of the recon- 
nection region and the magnetic field being absolutely 
straight, LV99 explores the real-world astrophysical re- 
connection. 

Our calculations in LV99 showed that the resulting 
reconnection rate is limited only by the width of the 
outflow region. This proposal, called "stochastic recon- 
nection", leads to reconnection speeds close to the tur- 
bulent velocity in the fluid. More precisely, assuming 
isotropically driven turbulence characterized by an in- 
jection scale, /, smaller than the current sheet length, we 
find 
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where m; is the velocity at the driving scale and Uturb is 
the velocity of the largest eddies of the strong turbulent 
cascade. Note, that here "strong" means only that the 
eddies decay through nonlinear interactions in an eddy 



turn over time (see more discussion of the LV99). All 
the motions are weak in the sense that the magnetic field 
lines are only weakly perturbed. 

It is useful to rewrite this in terms of the power injec- 
tion rate P. As the perturbations on the injection scale 
of turbulence are assumed to have velocities m; < Va, 
the turbulence is weak at large scales. Therefore, the re- 
lation between the power and the injection velocities are 
different from the usual Kolmogorov estimate, namely, 
in the case of the weak turbulence P 
Thus we get. 
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where / is the length of the turbulent eddies parallel to 
the large scale magnetic field lines as well as the injec- 
tion scale. 

The reconnection velocity given by equation (O does 
not depend on resistivity or plasma effects. Therefore 
for sufficiently high level of turbulence we expect both 
coUisionless and collisional fluids to reconnect at the 
same rate. 

2. Testing of Lazarian & Vishniac 99 Model 

Here we describe the results of a series of three di- 
mensional numerical simulations aimed at adding turbu- 
lence to the simplest reconnection scenario and testing 
equation We take two regions with strongly differ- 
ing magnetic fields lying next to one another The sim- 
ulations are periodic in the direction of the shared field 
(the z axis) and are open in the reversed direction (the x 
axis). The external gas pressure is uniform and the mag- 
netic fields at the top and bottom of the box are taken to 
be the specified external fields plus small perturbations 
to allow for outgoing waves. The grid size in the sim- 
ulations varied from 256x5 12x256 to 512x1028x512 so 
that the top and bottom of the box are far away from the 
current sheet and the region of driven turbulence around 
it. At the sides of the box where outflow is expected the 
derivatives of the dynamical variables are set to zero. 
A complete description of the numerical methodology 
can be found in Kowal et al. (2009). All our simula- 
tions are allowed to evolve for seven Alfven crossing 
times without turbulent forcing. During this time they 
develop the expected Sweet-Parker current sheet config- 
uration with slow reconnection. Subsequently we turn 
on isotropic turbulent forcing inside a volume centered 
in the midplane (in the xz plane) of the simulation box 
and extending outwards by a quarter of the box size. 
The turbulence reaches its full amplitude around eight 
crossing times and is stationary thereafter. 
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Figure 2: Left panel: Current intensity and magnetic field configuration during stochastic reconnection. We show a slice through the middle of the 
computational box in the xy plane after twelve dynamical times for a typical ran. The shared component of the field is perpendicular to the page. 
The intensity and direction of the magnetic field is represented by the length and direction of the arrows. The color bar gives the intensity of the 
current. The reversal in B, is confined to the vicinity of y=0 but the current sheet is strongly disordered with features that extend far from the zone 
of reversal. Right panel: Representation of the magnetic field in the reconnection zone with textures. 



In Figure 2 we see the current density on an xy slice 
of the computational box once the turbulence is well de- 
veloped. As expected, we see that the narrow stationary 
current sheet characteristic of Sweet-Parker reconnec- 
tion is replaced by a chaotic structure, with numerous 
narrow peaks in the current density. Clearly the pres- 
ence of turbulence has a dramatic impact on the struc- 
ture of the reconnection zone. In addition, we see nu- 
merous faint features indicating weak reconnection be- 
tween adjacent turbulent eddies. 

The speed of reconnection in three dimensions can 
be hard to define without explicit evaluation of the mag- 
netic field topology. However, in this simple case we 
can define it as the rate at which the x component of 
the magnetic field disappears. More precisely, we con- 
sider a yz slice of the simulation, passing through the 
center of the box. The rate of change of the area inte- 
gral of — Bx — is its flux across the boundaries of the 
box minus the rate at which flux is annihilated through 
reconnection (see more discussion in Kowal et al. 2009) 



where electric field is ^ - vx 1^ - rjj, Bx,exi is the abso- 
lute value of Bx far from the current sheet and L- is the 
width of the box in the z direction. This follows from the 
induction equation under the assumption that the turbu- 
lence is weak to lead to local field reversals and that the 
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Figure 3: Reconnection speed versus input power for the driven tur- 
bulence. We show the reconnection speed, defined by equation (4) 
plotted against the input power for an injection wavenumber equal to 
8 (i.e. a wavelength equal to one eighth of the box size) and a resis- 
tivity y„. The dashed line is a fit to the predicted dependence of P''^ 
(see eq. (3)). The horizontal line shows the laminar reconnection rates 
for each of the simulations before the turbulent forcing started. Here 
the uncertainty in the time averages are indicated by the size of the 
symbols and the variances are shown by the error bars. 
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Figure 4: Reconnection speed versus resistivity. We show the recon- 
nection speed plotted against the uniform resistivity of the simulation 
for an injection wavenumber of 8 and an injected power of one. We 
include both the laminar reconnection speeds, using the hollow sym- 
bols, fit to the expected dependence of ;;„, and the stochastic recon- 
nection speeds, using the filled symbols. As before the symbol sizes 
indicate the uncertainty in the average reconnection speeds and the 
eiTor bars indicate the variance. We included simulations with large, 
6j = 1, and small, B. = 0.1, guide fields. 

Stresses at the boundaries are weak to produce signifi- 
cant field bending there. In other words, fields in the x 
direction are advected through the top and bottom of the 
box, and disappear only through reconnection. Since 
we have assumed periodic boundary conditions in the 
z direction the boundary integral on the right hand side 
is only taken over the top and bottom of the box. By 
design this definition includes contributions to the re- 
connection speed from contracting loops, where Ohmic 
reconnection has occurred elsewhere in the box and 
decreases as the end of a reconnected loop is pulled 
through the plane of integration. It is worth noting that 
this estimate is roughly consistent with simply measur- 
ing the average influx of magnetic field lines through the 
top and bottom of the computational box and equating 
the mean inflow velocity with the reconnection speed. 
Following equation (|4| we can evaluate the reconnec- 
tion speed for varying strengths and scales of turbulence 
and varying resistivity. 

In Figure [3] we see the results for varying amounts 
of input power, for fixed resistivity and injection scale 
as well as for the case of no turbulence at all. The line 
drawn through the simulation points is for the predicted 
scaling with the square root of the input power The 
agreement between equation Q and Figure[3]is encour- 
aging but does not address the most important aspect of 
stochastic reconnection, i.e. its insensitivity to r]. 

In Figure |4] we plot the results for fixed input power 
and scale, while varying the background resistivity. In 



this case t] is taken to be uniform, except near the edges 
of the computational grid where it falls to zero over 
five grid points. This was done to eliminate edge ef- 
fects for large values of the resistivity. We see from the 
Figure |4] that the points for laminar reconnection scale 
as -\/?7, the expected scaling for Sweet-Parker reconnec- 
tion. In contrast, the points for reconnection in a turbu- 
lent medium do not depend on the resistivity at all. In 
summary, we have tested the model of stochastic recon- 
nection in a simple geometry meant to approximate the 
circumstances of generic magnetic reconnection in the 
universe. Our results are consistent with the mechanism 
described by LV99. The implication is that turbulent 
fluids in the universe including the interstellar medium, 
the convection zones of stars, and accretion disks, have 
reconnection speeds close to the local turbulent veloc- 
ity, regardless of the local value of resistivity. Magnetic 
fields in turbulent fluids can change their topology on a 
dynamical time scale. 

In Kowal et al. (2009) we also studied the depen- 
dence of the reconnection on the anomalous resistiv- 
ity, which increases effective resistivity for high current 
densities. The anomalous resistivity can be used as a 
proxy for plasma effects, e.g. collisionless eff'ects in re- 
connection. While it enhances the local speed of indi- 
vidual reconnection events results in Kowal et al. (2009) 
testify that the total reconnection rate does not change. 

Any numerical study has to address the issue of the 
possible numerical effects on the results. We show the 
dependence of the reconnection rate on the numerical 
resolution in Figure |5] The reconnection rate increases 
with the increase of the resolution, which testifies that 
the fast reconnection is not due to numerical effects. 
Indeed, higher numerical reconnection is expected for 
lower resolution simulations. 

Finally, it is important to give a few words in re- 
lation to our turbulence driving. We drive our turbu- 
lence solenoidally to minimize the effects of compres- 
sion, which does not play a role in LV99 model. The 
turbulence driven in the volume around the reconnec- 
tion layer corresponds to the case of astrophysical tur- 
bulence, which is also volume-driven. On the contrary, 
the case of the turbulence driven at the box bound- 
aries would produce spatially inhomogeneous imbal- 
anced turbulence for which we do not have analytical 
predictions (see discussion of such turbulence in Beres- 
nyak & Lazarian 2009). We stress, that it is not the 
shear size of our numerical simulations, but the corre- 
spondence of the observed scalings to those predicted 
in LV99 that allows us to claim that we proved that the 
3D reconnection is fast in the presence of turbulence. 
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Figure 5: Dependence of the reconnection rate on the numerical res- 
olution. If the fast reconnection were due to yet unclear numerical 
effects on small scales, we would expect to see the increase of the re- 
connection rate with the decrease of the numerical box. If anything, 
the actual dependence of the reconnection rate on the box size shows 
the opposite dependence. 
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Figure 6: Cosmic rays spiral about a reconnected magnetic field Hne 
and bounce back at points A and B. The reconnected regions move 
towards each other with the reconnection velocity Vr. The advec- 
tion of cosmic rays entrained on magnetic field lines happens at the 
outflow velocity, which is in most cases of the order of Va ■ Bounc- 
ing at points A and B happens because either of streaming instability 
induced by energetic particles or magnetic turbulence in the reconnec- 
tion region. In reality, the outflow region gets filled in by the oppo- 
sitely moving tubes of reconnected flux which collide only to repeat 
on a smaller scale the pattern of the larger scale reconnection. From 
Lazarian (2005). 



3. Acceleration of Cosmic Rays 

In what follows we discuss the first order Fermi accel- 
eration which arises from volume-filling reconnectior^. 



The LV99 presented such a model of reconnection and 
observations of the Solar magnetic field reconnection 
support the volume-filled idea (Ciaravella & Raymond 
2008). 

Figure |6] exemplifies the simplest realization of the 
acceleration within the reconnection region expected 
within LV99 model. As a particle bounces back and 
forth between converging magnetic fluxes, it gains en- 
ergy through the first order Fermi acceleration described 
in de Gouveia dal Pino & Lazarian (2003, 2005, hence- 
forth GL05) (see also Lazarian 2005). 

To derive the energy spectrum of particles one can 
use the routine way of dealing with the first order Fermi 
acceleration in shocks (see Longair 1992). Consider the 
process of acceleration of Mq particles with the initial 
energy Eq. If a particle gets energy ^Eq after a collision, 
its energy after m collisions is ff'Eo. At the same time if 
the probability of a particle to remain within the acceler- 
ating region is P, after m collisions the number of parti- 
cles gets P^Mo. Thus ln(M/M())/ ln(£'/£o) = Inf/luyS 
and 

Mo \Eq} 

For the stationary state of accelerated particles the num- 
ber M is the number of particles having energy equal or 
larger than E, as some of these particles are not lost and 
are accelerated further Therefore: 

const xE-'^''-''I"'PUE 



N(E)dE : 



(6) 



To determine P and /3 consider the following pro- 
cess. The particles from the upper reconnection re- 
gion see the lower reconnection region moving toward 
them with the velocity 2Vr (see Figure [3]). If a parti- 
cle from the upper region enters at an angle into the 
lower region the expected energy gain of the particle is 
6E/E - 2VrCos0/c. For isotropic distribution of par- 
ticles their probability function is p{0) = 2 sin cos 0d0 
and therefore the average energy gain per crossing of 
the reconnection region is 

r.n/2 



Vr 

c 



{6E/E}^— loos'- sin 0d0 ■■ 



4/3— (7) 

c 



An acceleration cycle is when the particles return back 
to the upper reconnection region. Being in the lower 
reconnection region the particles see the upper recon- 
nection region moving with the speed Vr. As a re- 
sult, the reconnection cycle provides the energy increase 
{6EIE),y,,, = %I3{VrIc) and 

P = ElEo = 1 + 8/3(y«/c) (8) 



'We would like to stress that Figure 1 exemplifies only the first 
moment of reconnection when the fluxes are just brought together. As 



the reconnection develops the volume of thickness A gets filled with 
the reconnected 3D flux ropes moving in the opposite directions. 
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Consider the case of Vrf,// <K Vr. The total number of 
particles crossing the boundaries of the upper and lower 
fluxes is 2 X 1 /4(nc), where n is the number density of 
particles. With our assumption that the particles are ad- 
vected out of the reconnection region with the magne- 
tized plasma outflow the loss of the energetic particles 
is 2 X Vrh. Therefore the fraction of energetic particles 
lost in a cycle is VRn/[l/4(nc)] - 4Vr/c and 

P = 1 - 4Vr/c. (9) 

Combining Eq. ([8]), (|9j one gets 

N{E)dE = constiE'^'^dE, (10) 

which is the spectrum of accelerated energetic particles 
for the case when the back-reaction is negligible (see 
also GLOSil 

The first order acceleration of particles entrained on 
the contracting magnetic loop can be understood from 
the Liouville theorem. As in the process of the mag- 
netic tubes are contracting, the regular increase of the 
particle's energies is expected. The requirement for the 
process to proceed efficiently is to keep the accelerated 
particles within the contracting magnetic loop. This 
introduces limitations on the particle diff'usivities per- 
pendicular to magnetic field direction. The subtlety of 
the point above is related to the fact that while in the 
first order Fermi acceleration in shocks magnetic com- 
pression is important, the acceleration via LV99 recon- 
nection process is applicable to incompressible fluids. 
Thus, unlike shocks, not the entire volume that shrinks 
for the acceleration, but only the volume of the magnetic 
flux tube. Thus high perpendicular diffusion of particles 
may decouple them from the magnetic field. Indeed, 
it is easy to see that while the particles within a mag- 
netic flux rope depicted in Figure 6 bounce back and 
forth between the converging mirrors and get acceler- 
ated, if these particles leave the flux rope fast, they may 
start bouncing between the magnetic fields of diff'erent 
flux ropes which may sometimes decrease their energy. 
Thus it is important that the particle diff'usion parallel 
and perpendicular magnetic field stays different. Parti- 
cle anisotropy which arises from particle preferentially 
getting acceleration in terms of the parallel momentum 
may also be important. 

4. Simulations of the Acceleration of Cosmic Rays 
by Reconnection 

In the numerical studies of the cosmic ray acceler- 
ation we use data cubes obtained from the models of 

'The obtained .spectral index is similar to the one of Galactic cos- 
mic rays. 



the weakly stochastic magnetic reconnection described 
in §2. For a given snapshot we obtain a full configura- 
tion of the plasma flow variables (density and velocity) 
and magnetic field. We inject test particles in such an 
environment and integrate their trajectories solving the 
motion equation for relativistic charged particles 

^^{jmu)^q[E + uxB), (11) 

where u is the particle velocity, y = (^l - u^/c^^ is the 
Lorentz factor, m and q are particle mass and electric 
charge, respectively, and c is the speed of light. 

The study of the magnetic reconnection is done using 
the magnetohydrodynamic fluid approximation, thus we 
do not specify the electric field E explicitly. Neverthe- 
less, the electric field is generated by either the flow of 
magnetized plasma or the resistivity effect and can be 
obtained from the Ohm's equation 

E^-vxB + 7]j, (12) 

where v is the plasma velocity and j = V x B is the 
current density. 

In our studies we are not interested in the acceleration 
by the electric field resulting from the resistivity effects, 
thus we neglect the last term. After incorporating the 
Ohm's law, the motion equation can be rewritten as 

j^iymit)^q[iit-xf)xB]. (13) 

In our simulation we do not include the particle en- 
ergy losses, so particle can gain or loose through the in- 
teraction with moving magnetized plasma only. For the 
sake of simplicity, we assume the speed of light 20 times 
larger than the Alfven speed Va, which defines plasma 
in the nonrelativistic regime, and the mean density is 
assumed to be 1 atomic mass unit per cubic centimeter, 
which is motivated by the interstellar medium density. 
We integrate equation [T3] for 10,000 particles with ran- 
domly chosen initial positions in the domain and direc- 
tion of the motion. 

In Figure|2]we show the particle energy evolution av- 
eraged over all integrated particles for two cases. In the 
left plot we used plasma fields topology obtained from 
the weakly stochastic magnetic reconnection models, in 
the right plot we use fields topology taken from the tur- 
bulence studies. Gray area shows the particle energy 
dispersion over the group of particles. 

In the case of reconnection model, the expected ex- 
ponential acceleration is observed until time about 100 
hours. Later on, the physical limitations of the compu- 
tational domain result in a different growth rate corre- 
sponding to £■ ~ f' "*^. In the case of turbulence without 
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Figure 7: Particle energy evolution averaged over 10,000 particles 
with the initial energy Eg = lO' MeV and random initial positions and 
directions. In the upper plot we show results for the weakly stochas- 
tic turbulence environment and in the lower plot for turbulent envi- 
ronment without magnetic reconnection. In both cases we assume 
c = 20Va and (p> = 1 u/cm^. 



Figure 8: Particle spectrum evolution for 10,000 particles with the 
uniform initial energy distribution £0 = 10^ - 10* MeV and random 
initial positions and directions. In the upper plot with show results for 
the weakly stochastic turbulence environment and in the lower plot for 
turbulent environment without magnetic reconnection. In both cases 
we assume c = 20Va and (p) = 1 u/cm^ . 



large scale magnetic reconnection the growth of energy 
is slower E ~ f ''^. This testifies that that the pres- 
ence of reconnection makes the acceleration more ef- 
ficient. The numerical confirmation of the first order 
acceleration in the reconnection regions is presented in 
our forthcoming paper. 

In Figure [8] we show the evolution of particle energy 
distribution. Initially uniform distribution of particle en- 
ergy ranging from 10^ to 10^ MeV evolves faster to 
higher energies if the reconnection is present. In this 
case, the final distribution is log-normal being more 
peaked over the time and with decreasing dispersion of 
energies in logarithmic scale. On the contrary, in the 
case of pure turbulence, the energy distribution after 
evolving to the log-normal shape preserves its disper- 
sion over the time. 

As magnetic reconnection is ubiquitous process, the 
particle acceleration within reconnection process should 
be widely spread. Therefore accepting the preliminary 
character of these results above we are involved in more 
extensive studies of the acceleration-via-reconnection. 



5. Explanation of the Anomalous Cosmic Ray Ori- 
gin 

The processes of the energetic particle acceleration in 
the process of turbulent reconnection can preaccelerate 
particles to the intermediate energies helping to solve 
the problem of particle injection into shocks. It can also 
act as the principal process of acceleration. Below we 
present the case where we believe that the latter takes 
place. 

Since the crossing of the termination shock (TS) by 
Voyager 1 (VI) in late 2004 and by Voyager 2 (V2) in 
mid 2007 it became clear that several paradigms needed 
to be revised. Among them was the acceleration of par- 
ticles. Prior to the encounter of the termination shock by 
VI the prevailing view was that anomalous cosmic rays 
(ACRs) were accelerated at the TS by diffusive shock 
acceleration (DSA) to energies 1-300 MeV/nuc (e.g., 
Jokipii & Giacalone, 1998; Cummings & Stone, 1998). 
However, with the crossing of the TS by VI the energy 
spectrum of ACR did not unroll to the expected source 
shape: a power-law at lower energies with a roll off at 
higher energies. After 2004, both the VI spectrum in 
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the heliosheath and the V2 spectrum upstream the TS, 
continued to evolve toward the expected source shape. 

To explain this paradox several models were pro- 
posed. Among them, McComas & Schwadron (2006) 
suggested that at a blunt shock the acceleration site for 
higher energy ACRs would be at the flanks of the TS, 
where the injection efficiency would be higher for DSA 
and connection times of the magnetic field lines to the 
shock would be longer, allowing acceleration to higher 
energies. Fisk et al. (2006) on the other hand sug- 
gested that stochastic acceleration in the turbulent he- 
liosheath would continue to accelerate ACRs and that 
the high-energy source region would thus be beyond the 
TS. Other works, such as Jokipii (2006) and Florinki 
and Zank (2006) try to explain the deficit of ACRs based 
on a dynamic termination shock. Jokipii (2006) pointed 
out that a shock in motion on time scales of the accel- 
eration time of the ACRs, days to months, would cause 
the spectrum to differ from the expected DSA shape. 
Florinki & Zank (2006) calculated the effect of Mag- 
netic Interacting Regions (MIRs) with the Termination 
Shock on the ACR spectral shape. They show that there 
is a prolonged period of depressed intensity in mid- 
energies from a single MIR. Other recent works have 
included stochastic acceleration, as well as other effects 
(Moraal et al., 2006, 2007; Zhang, 2006; Langner and 
Potgieter, 2006; Ferreira et al., 2007). It became clear 
after the crossing of the TS by V2 that these models 
would require adjustments. The observations by V2 
indicate for example that a transient did not cause the 
modulation shape of the V2 spectrum at the time of 
its TS crossing. When both spacecraft were in the he- 
liosheath in late 2007, the radial gradient in the 13-19 
MeV/nuc ions did not appear to be caused by a tran- 
sient. The 60-74 MeV/nuc ions have no gradient, so 
no north-south or longitudinal asymmetry is observed 
in the ACR intensities at the higher energies. 

In Lazarian & Opher (2009, LO09) we propose an 
alternative model, which explains the source of ACRs 
as being in the heliosheath and we appeal to magnetic 
reconnection as a process that can accelerate particles. 
LO09 explained the origin of the magnetic field rever- 
sals that induce magnetic reconnection in heliosheath 
and heliopause. 

Indeed, it is well known that magnetic field in the 
heliosphere change polarity and induce reconnection. 
For instance, as the Sun rotates magnetic field twists 
into a Parker spiral (Parker 1958) with magnetic fields 
separated by a current sheet (see Schatten 1971). The 
changes of magnetic field are also expected due to the 
Solar cycle activity. 

The question now is at what part of the heliosheath 



we expect to see reversals. The structure of the mag- 
netic field in the solar wind is complex. The solar mag- 
netic field lines near the termination shock are azimuthal 
and form a spiral (see Figure |9l). We expect the recon- 
nection and the corresponding energetic particle accel- 
eration to happen at the heliosheath closer to the he- 
liopause. This explains why Voyagers do not see the 
signatures of anomalous cosmic ray acceleration as they 
pass the termination shock. Appealing to their model of 
collisionless reconnection, Drake et al. (2010) provided 
a similar explanation of the origin of the anomalous cos- 
mic rays. 




Figure 9: Upper plot. Global view of the interaction of the solar wind 
with the interstellar wind. The spiral solar magnetic field (shown in 
dark dashed lines) is shown being deflected at the heliopause. The 
heliopause itself is being deflected by the interstellar magnetic field, 
(figure adapted from S. Suess (2006). Lower plot. A meridional view 
of the boundary sectors of the heliospheric cun'eiity sheet and how the 
opposite sectors get tighter closer to the heliopause. The thickness of 
the outflow regions in the reconnection region depends on the level of 
turbulence. From LO09. 
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6. Convergence with other Models of Reconnection 
and Acceleration 

Since the introduction of the LV99 model, more tra- 
ditional approaches to reconnection have been changed 
considerably. At the time of its introduction, the mod- 
els competing with LV99 were some modifications of a 
single X-point collisionless reconnection. Those mod- 
els had point-wise localized reconnection region and in- 
evitably prescribed opening of the reconnection region 
upon the scales comparable to L (see Figure 1). Such re- 
connection was difficult to realize in astrophysical con- 
ditions in the presence of random forcing which at high 
probability would collapse the opening of the reconnec- 
tion layer. Single X-point reconnection were rejected in 
observations of solar flares by Ciaravella & Raymond 
(2008). 

Modern models of collisionless reconnection resem- 
ble the original LV99 model in a number of respects. 
For instance, they discuss, similarly to the LV99, the 
volume filled reconnection, although one may still won- 
der how this volume filling is being achieved in the pres- 
ence of a single reconnection layer (see Drake et al. 
2006). While the authors still talk about islands pro- 
duced in the reconnection, in three dimensions these is- 
lands are expected to evolve into contracting 3D loops 
or ropes (Daughton et al. 2008), which is similar to what 
is depicted in Figure |6] Thus we do not expect to see a 
cardinal difference between the first order Fermi pro- 
cesses of the acceleration described in GL03 and later 
in Drake et al. (2006). This suggests that the backre- 
action of the particles calculated in Drake et al. (2006) 
considering the firehose instability may be employed as 
a part of the acceleration process described in GL03. 

The departure from the idea of regular reconnection 
and introduction of magnetic stochasticity is also obvi- 
ous in a number of the recent papers appealing to the 
tearing mode instabilit}0 as the process of enhancing 
reconnection (Loureiro et al. 2009, Bhattacharjee et al. 
2009). The 3D loops that should arise as a result of this 
process should be able to accelerate energetic particles 
via the process described in GL03. As tearing modes 
can happen in a coUisional fluid, this may potentially 



The idea of appealing to the tearing mode as a means of en- 
hancing the reconnection speed can be traced back to Strauss (1988), 
Waelbroeck (1989) and Shibata & Tanuma (2001). LV99 showed that 
the linear growth of tearing modes is insufficient to obtain fast recon- 
nection. The new attack on the problem assumes that the non-linear 
growth of the islands due to merging provides their growth rates at the 
large scales that are larger than the direct growth of the tearing modes 
at those scales. This situation when the non-linear growth is faster 
than the linear one is rather unusual and requires further investigation. 



open another channel of reconnection in such fluid. The 
limitation of this process is that the tearing mode re- 
connection should not be too fast as this would present 
problems with explaining the accumulation of the flux 
prior to the flare. At the same time the idea of tearing 
reconnection does not have the natural valve of enhanc- 
ing the reconnection speed, which is contrary to the LV 
model where the degree of reconnection is determined 
by the level of turbulence. Thus the periods of slow re- 
connection in LV99 model are ensured by the low level 
of turbulence prior to the flare. We believe that tearing 
reconnection can act to distabilize the Sweet-Parker re- 
connection layer, inducing turbulence. 

We note, however, that in most astrophysical situa- 
tions one has to deal with the pre-existing turbulence, 
which is the consequence of high Reynolds number of 
the fluid. Such turbulence may modify or suppress in- 
stabihties, including the tearing mode instability. We 
claim that it, by itself, induces fast reconnection. We 
may also argue that even if the astrophysical fluid is 
kept initially laminar the fluid a thick outflow from the 
reconnection region caused by tearing is expected to be- 
come turbulent. We suspect that this may be the cause of 
the reconnection explosions reported recently (Lapenta 
2008, Bettarini & Lapenta 2009). 

7. Summary 

The successful testing of the LV99 model of fast re- 
connection opens avenues for the search of implications 
of that scheme. One of the implications of the model is 
the first order Fermi acceleration of energetic particles 
in the reconnection layer As reconnection processes 
are expected to be ubiquitous in astrophysics, we expect 
the acceleration in reconnection layers to be also ubiqui- 
tous. Our simulations of the energetic particle accelera- 
tion in the reconnection layer provide results consistent 
with the first order Fermi acceleration. The origin of the 
anomalous cosmic rays may be related with the mecha- 
nism of particle acceleration via reconnection. 
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